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Charge Modulations vs. Strain Waves in Resonant X-Ray Scattering 
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A method is described for using resonant x-ray scattering to separately quantify the charge (va- 
lence) modulation and the strain wave associated with a charge density wave. The essence of the 
method is a separation of the atomic form factor into a "raw" amplitude, fn{uj), and a valence- 
dependent amplitude, fo(cj), which in many cases may be determined independently from absorp- 
tion measurements. The advantage of this separation is that the strain wave follows the average 
quantity |/_r(w)+ <v> fr>(u>)\ 2 whereas the charge modulation follows only \fn(u)\ 2 - This allows 
the two distinct modulations to be quantified separately. A scheme for characterizing a given CDW 
as Peierls-like or Wigner-like follows naturally. The method is illustrated for an idealized model of 
a one-dimensional chain. 

PACS numbers: 78.70.Ck, 87.64.Gb 
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Charge density waves (CDWs) are pervasive in con- 
densed matter, arising in di- and trichalcogenides, con- 
ducting polymers, transition metal oxides etc. However 
they may form by a variety of mechanisms. The best 
understood is the Peierls mechanism in which a gap at 
the Fermi level opens, and a charge (valence) modulation 
forms, through the presence of a strain wave. The basis 
for the Peierls mechanism is the clectron-phonon interac- 
tion which drives pairing between an electron and a hole 
with equal but opposite wave vector, kp, resulting in a 
charge modulation with total momentum 2kp. A CDW 
may also be driven by electron interactions, however, an 
idealized example being the Wigner crystal (WC) ground 
state of a dilute electron gas, which is driven by di- 
rect Coulomb0. Because of the close analogy between 
charge density waves and superconductivity |2|, and the 
occurance of both in many systems such as the dichalco- 
genides and the copper-oxides in which superconductivity 
is not clearly understood, it would be extremely useful 
to distinguish between Peierls-like and Wigner-like (i.e. 
interaction-driven) mechanisms in practice[3|- 

The main distinction between a Peierls CDW and a 
WC is the presence of a lattice distortion in the former 
case that is not required in the latter. However in any 
real material even a WC must have at least some inci- 
dental lattice distortion if only because of electrostatics. 
Therefore, in practice, distinguishing between these two 
phenomena is a quantitative rather than a qualitative 
matter. 

X-ray diffraction is sensitive to both charge modula- 
tions and strain waves. Scattering from both can be en- 
hanced, through different resonant mechanisms, by tun- 
ing the x-ray photon energy to an atomic core transition, 
i.e. so-called resonant x-ray scattering (RXS). Here a 
general method is proposed for using RXS to distinguish 
between scattering from a charge modulation and scatter- 
ing from a strain wave, even when the two coexist. This 
method allows the charge and strain amplitudes to be 
quantified independently. A system classifying a given 
charge density wave as either Peierls or Wigner in ori- 
gin naturally arises from this separation. The method is 



demonstrated for the idealized case of a one-dimensional 
chain. 



I. SEPARATING THE ATOMIC FORM FACTOR 

In quantum electrodynamics, the cross section for elas- 
tic scattering of photons from a material ispj 
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where p q is a Fourier component of the total electron 
density, q is the momentum transfer, p is the momen- 
tum operator, A is the photon field operator (vector po- 
tential), and u> is the photon energy. The first of these 
two terms is equivalent to classical Thomson scatter- 
ing, does not depend on u, and gives rise to the "nor- 
mal" dispersion of optical constants in matter, i.e. in SI 
units x q = — r e A 2 p q /47T 2 , where r e is the classical elec- 
tron radius and A = 2ttc/lu 5]. The second, resonant 
term corresponds to scattering via virtual transitions be- 
tween core and valence states, is highly dependent on 
u>, and is the origin of anomalous dispersion 6] . This 
term, sometimes referred to as the Kramers-Heisenberg 
formula, has been shown to provide direct sensitivity to 
valence band charge and spin ordering in condensed mat- 
ter |Sllll|l3|ll|ll|ll|ll|lllillI3, provided w 
coincides with an edge threshold. 

How can one distinguish between a charge moduation 
and a strain wave? One might be tempted to just com- 
pare the relative sizes of nonresonant and resonant scat- 
tering on the assumption that a lattice distortion, which 
translates all the electrons in an atom including those in 
the core, contributes mainly to the former and valence 
band effects mainly to the latter. Indeed, the Kramers- 
Heisenburg formula alone is normall y u sed to analyze 
near-edge scattering data [lfj, IllL Il4| . However one 
must realize two things. First, the valence modulation it- 
self should scatter in the Thomson channel as well since 
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it corresponds to a real, albeit small, charge modulation. 
Second, it was realized long ago 18J that the resonant 
term also contributes to scattering far from the edges via 
virtual, off-shell processes. The optical contants of Cu 
metal, for example, exhibit no region of 'normal' disper- 
sion from the infrared up to 100 keV[l8l|. So the resonant 
term is omnipresent, a trait that is commonly exploited in 
crystallographic phasing methods[19|. More to the point, 
both the valence and lattice modulations of a CDW scat- 
ter in both channels, and a more appropriate division of 
the cross section than Eq. 1 is needed. 

Here it is proposed to divide the cross section instead 
into a 'raw' and 'valence-dependent' component. More 
specifically, assuming the form factor for an atom, /, is a 
function of both u> and the valence state, v, it is proposed 
to make the separation 

f(w,v) = f R (u)+vf D (w). (2) 

Here Jr is the 'raw' part of the atomic scattering factor 
and fr> describes how it changes with valence 20]. Like 
the full /, Jr is dimensionless (i.e. 'electrons') and con- 
verges to Z* as ui — > oo@- Jr contains both resonant 
and nonresonant processes but only those, such as edge 
jumps, that are independent of the detailed electronic 
structure of the atom. 

Jd has units of electrons/valence and can be thought 
of the scattering power of the valence modulation. Like 
Jr this quantity contains both resonant and nonresonant 
processes, however the nonresonant contribution will be 
extremely small. One can expect fo to be large only for 
to near the threshold of an absorption edge, where the 
intermediate states involve transitions directly into the 
valence band. 

There are two advantages to the division in Eq. 2. 
First, as will be illustrated in the next section, scattering 
from a strain wave will follow only <w>)| 2 , where 
<u> is the average atomic valence of the material, while 
scattering from the valence modulation follows |/d(u;)| 2 . 
By measuring the photon energy dependence of a CDW 
reflection, provided ju and ]r are known, one can deter- 
mine separately the charge and strain amplitudes as well 
as their relative phase. 

Second, and most importantly, in many cases f R 
and fn can be independently determined from doping- 
dependent x-ray absorption (XAS) data. More specifi- 
cally, / is related to the absorption coefficient [j| by 



Im[n(uj, <«>)] = — ^r-ri-Im [f R (u>)+ <v> f D (uj)} (3) 
ZirVc 

where n is the complex index of refraction, N is the 
number of atoms in a unit cell, and V c is the unit cell 
volume. So, by measuring the absorption coefficient on 
two or more samples of different <v>, Im\f R {u})\ and 
Im[fr)(uj)] can be determined by solving a system of 
equations at each value of u. The real parts can then 
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FIG. 1: One-dimensional, monatomic chain. u n and v n de- 
scribe respectively the displacement from equilibrium and va- 
lence state of the nth atom in the lattice. 



by determined by Kramers-Kronig transform. Without 
having to appeal to a specific model, then, the charge and 
strain amplitudes can be determined. XAS acts as a ref- 
erence against which resonant scattering measurements 
can be calibrated. 

Use of this procedure rests on several assumptions. 
First, one must assume that the scattering processes are 
local and that an atomic form factor, /, is definable. This 
will be valid as long as v g r <C A/ 2, where v g is the group 
velocity of the core electron-hole pair, r is its radiative 
lifetime, and A is the x-ray wavelength. For measure- 
ments near edges with sharp white lines, such as the L 
edges of the transition metals, one can expect this con- 
dition to be well satisfied. Next one must assume that 
changes in XAS spectra due to changes in <i>>, de- 
termined by comparing samples with different chemical 
composition, indeed arise from valence effects and not 
extrinsic phenomena such as changing crystal structures 
etc. Next, one must assume that changes in / are small 
i.e. are linear in v, which should be true as long as the 
CDW amplitude is not too large. It is not necessary that 
the sample have a rigid or noninteracting band struc- 
ture. In determining f R and fo from XAS one implic- 
itly makes the assumption that the electronic structure 
of every point in the CDW corresponds to the average 
electronic structure in a sample with that valence. Fi- 
nally, one assumes that x-ray absorption measurements 
are a good measure of the forward scattering amplitude, 
i.e. that there are not hidden loss processes such as pho- 
toelectron production. Violation of this last point has 
already been observed but is not a significant effect. 



II. ONE-DIMENSIONAL CHAIN 

To see how f R and fo enter a resonant scattering ex- 
periment consider the idealized case of a one-dimensional, 
monatomic chain, as shown in Fig. 1. Imposed upon this 
chain are a longitudinal strain wave, described by a set 
of displacements u n , and a valence modulation, described 
by a set of valence states, v n . For simplicity it is assumed 
that the modulations are sinusoidal, i.e. 



U n — uo cos(fcna) (4) 



3 



v n = <v> + vo cos(kna + 4>) (5) 

where u n implies that the nth atom resides at position 
r n = na+u n , where a is the average lattice parameter, uq 
and vq are the aplitudes of the strain and charge waves, 
respectively, <i>> is the average atomic valence, and k is 
the wave vector of the CDW. The quantity <f> allows for 
the fact that, while the two modulations have the same 
wave vector, there may be a phase difference between 
them. 

The integrated intensity of a Bragg reflection is pro- 
portional to the square of the scattering amplitude, p q , 
which is given by 



^) = ^£/(<^»)e^" (6) 
n=l 

where f(u>,v n ) is as defined by eq. (2) and N — 2ir/ka is 
the number of sites in the supercell, i.e. the dimensionless 
CDW period [2l| • q is the momentum transfer along the 
chain and is restricted to discrete values q = 2irl/Na, 
where I is an integer. If one multiplies out the terms one 
arrives at four distinct contributions to the scattering 
amplitude. 

Pq {u) = p B q {uj) + p\{u>) + />» + p™(u) (7) 
The first term, 



p?M = -/(*,<«>) 5>^, (8) 

71=1 

is the "Bragg" term and corresponds to resonant x- 
ray scattering off the average lattice. This quantity 
is independent of uq and vq and is nonzero whenever 
q = 2nm/a, where m is an integer, i.e. at the Bragg 
points of the undistorted chain. This term always has 
exactly the value f(uj,<v>), which demonstrates that 
regular Bragg scattering has the energy dependence of 
\f(ui, <v>)\ 2 = \f R (cu)+ <v> f D (uj)\ 2 , i.e. simply tracks 
the average scattering factor of the atomic lattice, pf is 
highly resonant but nonzero for all values of u). 
The next term 



= - v f D (w) £ cos(fcna + fie**™ (9) 

n— 1 

is the "valence" term and corresponds to resonant x- 
ray scattering off the valence modulation. This quan- 
tity is nonzero only for q = ±fc and has the value 
p v k (w) = Uq/d(w) exp(i</>)/2. This demonstrates that res- 
onant scattering from the valence modulation tracks only 
|/d(w) 1 2 and is proportional to Vq. Since fu is significant 



only near an edge, one can expect to see scattering from 
a valence modulation only over a narrow energy range 
near threshold. This is the experimental signature of a 
valence modulation in resonant x-ray scattering. 

The third term is the "strain" term and corresponds 
to scattering off the lattice distortion. Provided the size 
of the distortion is small, i.e. cxp(iqr n ) = exp(iqna)(l + 
iqu n ), the strain term has the value 
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P» = j^iquof(cJ,<v>) £cos(fcna)e l «" a . (10) 

n=l 

This quantity is also nonzero only if q = ±fc and reduces 
to Pk(u>) = iquof(u>, <v>)/2. Scattering from the strain 
wave is proportional to (fcuo) 2 , as expected, and like the 
Bragg term follows \f(u>, <u>)| 2 so is visible at all values 
of uj. Evidently all structural scattering is alike in its 
adherence to the average f(u>,<v>). Unlike the Bragg 
term, however, the strain term occurs at the same q as 
the valence scattering and and the two may coherently 
interfere. We will see that, if this interference is visible, 
it provides a means to determine the phase shift, <f>. 

The final term is somewhat unexpected and apparently 
has not been addressed before. It is a mixed term corre- 
sponding to coherent displacement of the valence modu- 
lation. Written out, it has the form 



Pq v (w) = —iqv u f D (Lj)^2[cos(2kna+4))+cos((t))}e lqr " 

n—1 

(11) 

This quantity is nonzero only if q = ±2fc and, while it 
arises from both strain and valence scattering, it tracks 
only fr>. p uv is not a multiple-scattering effect; it is a sign 
that if both charge and strain modulations are present 
the total modulation is anharmonic. Observation of this 
term, i.e. by tuning the x-ray energy near threshold and 
scanning around 2k, would be a strong validation of our 
approach. However this term is extremely small. 



III. A GENERIC EDGE 

The various scattering processes are best illustrated 
with a specific model of /r and fn near an absorption 
edge for a single atom. As was argued earlier, fn de- 
scribes processes that are material- and valence- indepen- 
dent, such as the edge jump, and fo describes only those 
processes that depend on the atomic valence, v. As an 
illustrative model of we consider a generic edge jump, 
i.e. a dielectric susceptibility whose imaginary part has 
the form 



Im[x(u)] = -9(u-u e )- (12) 
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FIG. 2: Real and imaginary parts fn and fn as determined 
from the parameter values chosen (see text). 



where uj e is the edge energy and J is the size of the jump. 
X must satisfy the Kramers-Kronig relations so its real 
part is given by the integral 
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(13) 

This can be considered a generic analytic model for an 
idealized edge. The atomic scattering factor is related to 
this susceptibility bv[22j 
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(14) 



where x is m C S S units. V c here is the volume per atom. 
Bear in mind that / in general is a tensor, but for sim- 
plicity we will treat it as a scalar here. It is important to 
point out that this quantity actually diverges as uj — » 00 
but near the edge it is well behaved. 

Eq. 14 is for an isolated edge. However in practice 
there are other transitions that are far away but still 
contribute to /. To account for these we add a constant 
"background" scattering factor, /o, i.e. 



/flM -> fn(uj) + fa 



(15) 



Using the notation f(u>) = f l (u>) + if 2 (uj), writing it out 
explicitly, fn has the form 
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-A log 
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f R (uj)=A—6(uj-uj e ) + f 2 

0J e 



(16) 



(17) 



where A = V c J/r e c 2 . 

The valence-dependent form factor, fjj, could take 
on many forms. In high temperature superconductors, 
for example, it actually exhibits a sign change [l7j. For 
present, illustrative purposes, we will simply take it to 
be a Lorentzian, i.e. 



fo(oj) 



A 



w - «7 



(18) 



In principle fa should have an energy-independent com- 
ponent representing Thomson scattering from the valence 
modulation, but in situations of interest this is small. 

For numerics, we will take the parameter values f® = 
9.50, / 2 ° = 0.246, A = 4.114, A = 2 eV, and 7 = 0.2 
eV, uje = 534 eV, and uj = 528.6 eV. While it is in- 
tended that this illustration be general these parameters 
are quite a good model of the oxygen K edge, fn and fr> 
for these parameter values are plotted in Fig. 2. 

For the CDW itself we use the parameters vq = 0.1 
electron, uq — 0.1a, k = 2ir/4a, and <t;>= 0.12. The 
scattered intensity at q = k is given by square of the total 
scattering amplitude, \p%(uj) + p v k {uj)\ 2 . The strain and 
charge scattering amplitudes can coherently interfere, so 
must be added before squaring. The resulting quantity 
depends explicitly on the phase, (j). 

The quantity \p^(uj) + p^(uj)\ 2 is plotted against energy 
in Fig. 3. For this plot cj> is taken to be zero. Spectra 
are shown for a pure charge wave, a pure strain wave, 
and a composite wave. Notice that scattering from the 
charge-only wave is nonzero only in the region in which 
/d, displayed in Fig. 2, is nonzero. The strain-only wave, 
in contrast, is visible at (almost) all energies. This affirms 
ones intuition that scattering from structural distortions 
should be visible at all energies, but the valence modula- 
tion only near the edge. If both modulations are present 
the lineshape is not simply the sum of the two because 
of the non-trivial dependence on the phase factor, <j>. 

To illustrate this phase dependence the line shape is 
plotted in Fig. 4 for vaious values of <\>. Notice that 
not only the intensity of the various features but in fact 
the entire spectral lineshape depends sensitively on 4>. 
Therefore, if both strain and charge scattering are si- 
multaneously visible, and fn and fn are determined in- 
dependently from XAS, it should be possible to objec- 
tively determine the phase from a one-parameter fit to 
this shape. 



IV. PEIERLS VS. WIGNER CDWS 

Distinguishing between a Peierls CDW and a more ex- 
otic CDW driven by many-body interactions, such as a 
stripe phase or a Wigner crystal, is a quantitative rather 
than a qualitative matter. This is because in any real 
material even a CDW driven purely by electron-electron 
interactions, because of electrostatics, must still be ac- 
companied by a lattice distortion, though it may be 
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FIG. 3: Plots of the integrated intensity of resonant x-ray 
scattering at q = k, i.e. \p^(uj) + pl(w)\ 2 . Here = 0. 
(black line) Charge modulation only, (red line) Strain wave 
only, (blue line) Both charge and lattice modulations present. 
Charge scattering is localized near threshold but the strain 
wave, while resonant, is visible at all energies. Notice that 
the lattice distortion has a strong intensity max at the edge 
jump due to the cusp in the real part of fn. 



small. Identifying a given CDW as either Peierls-like or 
Wigner-like requires a quantitative comparison between 
the charge and lattice amplitudes. Once uo and vo have 
been determined by the procedure just outlined, then, a 
quantity of great interest is the ratio 



(19) 



This quantity has units of length -1 and describes the 
degree to which the CDW is Wigner-like, i.e. diven by 
many-body interactions rather than a Peierls distortion. 
For example, a perfect Wigner crystal with no lattice 
distortion at all would have W — oo. A typical Peierls 
CDW on the other hand, such as that in NbSe3 , has a lat- 
tice distortion of approximately uq ~ O.OlA and a charge 
amplitude of vq ~ 0.1, giving W = 10A -1 . We propose 
that a CDW with W less than about 20 should be con- 
sidered a Peierls CDW. If W > 100 the CDW probably 
arises at least partly from many-body effects and should 
be considered "exotic" . It is likely that CDWs can exist 
over the entire continuum of values of W . It would be 
particularly enlightening to determine the W values for 
several CDW materials, such as the copper-oxides and 
the dichalcogenides, that also exhibit superconductivity. 
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FIG. 4: Plots of the integrated intensity of resonant x-ray 
scattering at q = k, i.e. \p^(ui) + p£(cj)| 2 for various values of 
the phase between charge and strain modulations, <f>. Notice 
that the spectral shape changes with (f>. If fo and Jr are 
known from absorption data, <j) can he determined through 
analysis of this shape. 



V. SUMMARY 

A descripton of resonant x-ray scattering was intro- 
duced in which the atomic scattering factor is divided 
into raw- and valence-dependent amplitudes, fn and fjj. 
The advantage of this division is that resonant x-ray scat- 
tering from the strain wave component of a CDW tracks 
the average form factor |/r(w)+ <t;> fo{u)\ 2 whereas 
the charge (valence) scattering tracks only \fu(uj)\ 2 . In 
many cases fn and /jj can be independently determined 
from x-ray absorption measurements on materials with 
different average valence, combined with Kramers-Kronig 
analysis. This provides a means to separately quantify 
the charge and strain components of a CDW. In this 
framework one can define a quantity "W" which pro- 
vides a quantitative means to characterize a given CDW 
as Peierls-like, Wigner-like, or anywhere on the contin- 
uum between. 
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